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We present a measurement of the longitudinal spin cross section asymmetry for deep inelastic 
muon-nucleon interactions with two high transverse momentum hadrons in the final state. Two 
methods of event classification are used to increase the contribution of the Photon Gluon Fusion 
process to above 30%. The most effective one, based on a neural network approach, provides the 
asymmetries A e p N ^ ehhx = 0.030 ± 0.057 ±0.010 and Af^ hhx = 0.070 ±0.076 ±0.010. From these 
values we derive an averaged gluon polarization AG/G = — 0.20±0.28±0.10 at an average fraction 
of nucleon momentum carried by gluons (rj) = 0.07. 

PACS numbers: 13.60.Hb, 13.88.+e, 14.70.Dj 



I. INTRODUCTION 

The Spin Muon Collaboration (SMC) has extensively 
studied polarized deep inelastic lepton-nucleon scattering 
using the high energy muon beam at CERN and large po- 
larized hydrogen and deuterium targets. This program 
was initiated by the observation in a previous CERN ex- 
periment (EMC) that only a small fraction of the proton 
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spin is carried by the spin of the quarks y|. The SMC 
results have confirmed this observation for protons and 
provided the first measurement of the spin structure of 
deuterons which allowed for the verification of the fun- 
damental Bjorken sum rule 

The high energy polarized data from SMC, combined 
with the high precision data from the SLAC and 
DESY Q experiments at lower energy, cover a kinematic 
range allowing for a QCD analysis of the spin struc- 
ture function g\. Various analyses have been performed 
at next-to-leading order with different input parameteri- 
zations for the polarized parton densities and different 
choices of the fitted parameters 0, Q. They provide 
consistent results for the polarized quark densities but 
bring little information on the polarized gluon density 
AG. This is an expected feature since g\ is sensitive to 
gluons only through its Q 2 evolution and the available g\ 
data cover only a narrow range in Q 2 at a given value of 
x. In particular it is still not possible to test the hypoth- 
esis, formulated many years ago, that the gluon spin may 
account for a sizable fraction of the nucleon spin . 

A direct measurement of the gluon polarisation is pos- 
sible via the Photon Gluon Fusion (PGF) process, which 
is illustrated in FigQ] together with the two other lowest 
order diagrams: the virtual photon absorption (leading 
process "LP" ) and gluon radiation (QCD Compton scat- 
tering "QCD-C"). Since the contribution of the PGF 
diagram is small, the event selection procedure should 
be very effective in discriminating the PGF process from 
other channels. This can be achieved either by select- 
ing events where a charmed particle is produced (e.g. a 
D meson) or events with hadrons of large transverse mo- 
menta (pt) relative to the virtual photon direction |9lllfj|. 
Both possibilities will be used in the COMPASS experi- 
ment presently running at CERN [llj. 

In this paper we present an evaluation of the gluon 
polarization, AG/G, from the SMC data. We limit the 
analysis to the DIS region (Q 2 > 1 GeV 2 ) and select 
events with high px hadrons. The SMC experimental 
setup was not optimized for the detection of hadrons 
produced at large angles, so the precision of the result 
is obviously limited. This is, however, the first attempt 
to tag PGF with light quark production in a DIS exper- 
iment. 

A determination of the gluon polarization from events 
with high pt hadrons has been attempted on the ep data 
from the HERMES experiment ^2 a t lower incident en- 
ergy and in a kinematic range where quasi-real photo- 
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production is dominant. 



II. FORMALISM 

Experimentally observed spin-dependent effects are 
small and have to be determined from the cross section 
asymmetry defined as the ratio of polarized (Act) and 
unpolarized (ct) cross sections 
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where ][ and ]] refer to anti-parallel and parallel con- 
figurations of the nucleon and incoming lepton spins. At 
the parton level the hard scattering cross section consists 
of three terms corresponding to the LP, QCD-C and PGF 
processes. According to the factorization theorem ct and 
Act can be written as convolutions of the parton distri- 
butions (F, AF), hard-scattering cross sections (ct, Act) 
and fragmentation functions D of partons into hadrons: 



ct = F ® a ® D 



Act = AF(g> Act® D. 



(2) 



The parton distributions stand for quarks, antiquarks, 
and gluons. The spin dependent distributions are de- 
noted Aq = q* — q+ for quarks, antiquarks and AG = 
G^ — G^ for gluons and the corresponding spin-averaged 
ones q = q 1 + and G = G^ + . Here, the up and 
down arrows correspond to parallel and anti-parallel con- 
figurations of the parton and nucleon spins. 

After insertion of the full expression for ct and Act into 
Eq. ft}, the final expression for the cross section asym- 
metry with production of at least two hadrons with large 
transverse momenta, A eN ^ ehhx , reads 



j^tN->ehhX 
Aq 
9. 

AG 
~G 



— ({d L L) LtJ RLP + {a L L) QCD - C RQCD-c) + 



{aLL) P F RpGF, 



(3) 



where = (Aa/2a) are the average partonic asym- 

metries and R the cross section ratios of the different pro- 
cesses shown in Fig. 2] with respect to the total cross sec- 
tion in the selected sample. The asymmetry A lN ^ e - hhx 
thus permits an evaluation of the gluon polarization if all 
other elements in Eq. |[3J| are known. The quark asymme- 
try Aq/q is approximated by the value of A\ obtained in 
inclusive measurements. The partonic asymmetries &ll 
are calculated for simulated events and averaged over the 
selected sample; in the kinematic region covered by the 
SMC data, they are found to be positive for the first two 
processes and negative for PGF. The ratios R are taken 
from the simulated sample to which the same selection 
criteria are applied as to the data. 



The statistical precision of the gluon polarization de- 
termined from Eq. ||2J) depends on the precision of the 
measured asymmetry A iN ^ lhhX anc j on ^ e fraction of 
PGF events (Rpgf) in the final sample. Therefore the 
aim of the present analysis is to select a large enough 
sample with a maximal contribution of PGF events. 

The description of hadron production in DIS muon 
data in terms of the three processes of Fig. [T] has 
been successfully tested in previous experiments [w| . 
Other processes, such as those involving resolved pho- 
tons, are expected to have small contributions for Q 2 
above 1 GeV 2 and are not considered here. 



III. THE EXPERIMENT 

The experimental setup at the CERN muon beam con- 
sisted of three major components: a polarized target, a 
magnetic spectrometer and a muon beam polarimeter. A 
detailed description of the experiment and of the analysis 
of the inclusive data can be found in Refs. H0. The 
muon beam polarization, Pg, was determined from the 
spin asymmetries measured in polarized muon-electron 
scattering and from the energy spectrum of positrons 
from muon decays and was found to be —0.795 ± 0.019 
for an average beam energy of 187.4 GeV [3. The 
target consisted of two cells filled with butanol, deuter- 
ated butanol or ammonia [l7| . The two cells were po- 
larized in opposite directions by dynamic nuclear polar- 
ization. The average target polarizations, Pt, were ap- 
proximately 0.90 for protons and 0.50 for deuterons, with 
a relative error APt/Pt of 3-5%. The polarization was 
reversed five times a day. 

The counting rate asymmetry, A exp , is determined 
from the number of events counted in upstream and 
downstream target cells before and after polarization re- 
versal. This is done by solving the resulting second order 
equation, as described in [l8|. 

The cross-section asymmetry, A eN ^ ehhx j j s related to 

A exp by . 



A 



IN^lhhX 



Pe.Pt/ 



Aexp i 



(4) 



where / is the effective dilution factor, which takes into 
account the dilution of spin asymmetries by the presence 
of unpolarizable nuclei in the target and also by radiative 
effects on the nucleon. The effect of unpolarizable materi- 
als can be expressed in terms of the numbers ua of nuclei 
with mass number A and the corresponding total spin- 
independent cross sections ct^ 3 '. The radiative effects 
on the nucleon p^L are taken into account through 
the ratio of one photon exchange to total cross-sections 
p = Vp J d /&p°d- The evaluation of the effective dilution 
factor for inclusive events and for events with observed 
hadrons is described in Ref. 0. Polarized radiative cor- 
rections are applied to the asymmetries as described in 
Refs. |2(j. In this analysis polarized radiative cor- 
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rections and dilution due to radiative effects are reduced 
because processes without hadrons are excluded. 



IV. SAMPLE SELECTION 

The total sample of data collected by the SMC ex- 
periment during the years 1993-1996 with muon beam 
of _E=190 GcV and longitudinally polarized target was 
used for the analysis. It consists of samples of similar 
size taken on polarized protons and deuterons. 

The standard cuts on inclusive kinematic variables 3] , 
v = E — E' >15 GeV and E >19 GeV were imposed to 
reject events with poor kinematic resolution and muons 
from hadron decay, respectively. The cut y = v/E <0.9 
removes a region where the uncertainty due to radiative 
corrections becomes large. Two other cuts were applied 
in close relation to the formalism used in the analysis: a 
cut Q 2 >1 GeV 2 rejects the region dominated by non- 
perturbative effects and allows to interpret the results in 
terms of partons. A cut y >0.4 removes events which 
carry little spin information due to a small virtual pho- 
ton polarization. In addition, cuts on the muon scat- 
tering angle were applied in order to match the angular 
acceptance of the hardware triggers. 

In the leading process (LP) most hadrons have small 
Pt as only the intrinsic kr of quarks in the nucleon pl| 
and the fragmentation mechanism contribute to it. A 
different situation occurs for QCD-C and PGF, where 
hadrons mainly acquire transverse momentum from pri- 
marily produced partons. For this reason, the require- 
ment of two observed hadrons with large transverse mo- 
menta enhances the contribution of the PGF and QCD-C 
processes in the selected sample. 

In the present analysis, the events of interest include 
a reconstructed beam muon, a scattered muon, and at 
least two charged hadrons. They represent about 20% 
of the total number of events with reconstructed beam 
and scattered muons, used for inclusive studies. Hadron 
tracks were accepted if they could be associated to the 
primary interaction point, i.e. the vertex, defined by the 
beam and scattered muons. The same association criteria 
as in the SMC analysis of Ref. were applied. In order 
to suppress the contribution from the target fragmenta- 
tion region, cuts on the reduced longitudinal momentum 
of the hadron, xf > 0.1, and on the hadron fractional 
energy, z — Ey t jv > 0.1, were applied. 

The further requirement of two hadrons with pr > 
0.7 GeV selects about 5% of the events passing all pre- 
vious cuts. The electron contamination to this sample is 
expected to be negligible because electrons are generally 
produced at low px- This is confirmed by the ratio of 
the energy deposited in the electromagnetic part of the 
calorimeter to the total deposited energy, which does not 
show any peak at 1.0 for tracks with pr > 0.5 GeV. Af- 
ter all selections the total number of remaining events 
amounts to about 80k for the proton and 70k for the 
deuteron sample. 



V. MONTE CARLO SIMULATION 
A. Conditions for MC generation 

The interactions were simulated using the LEPTO 6.5 
generator |22| with a leading order parameterization of 
the unpolarizcd parton distributions 23]. The spin de- 
pendent effects were calculated using POLDIS 24| with a 
consistent set of polarized parton distributions |25| . The 
kinematic limits of the MC generation were defined so as 
to cover the full kinematic region of the data. Default 
values were used for most of the steering parameters of 
the LEPTO generator. Below we discuss only the modi- 
fied conditions and parameters. 

The matrix elements of first order QCD processes ex- 
hibit collinear divergences in the cross channel and dif- 
ferent schemes are used to avoid such singularities. The 
so-called zs scheme, which allows for lower values of the 
7*-parton center of mass energy y§, was used in the 
simulation with modified cut-off parameters. The effect 
of the cut-off values on any observable distribution for 
events with high px hadrons is only marginal. 

The description of interactions requires the choice of 
two scales: a factorization scale, which appears in the 
parton densities, and a renormalization scale which ap- 
pears in expressions depending on the strong coupling 
constant a s . Here the usual choice of Q 2 was made in 
both cases. In these conditions, after kinematic cuts on 
event variables only, the generated sample contains 8% 
PGF events. 

In order to describe the data, it was found necessary 
to change the values of two fragmentation parameters in 
JETSET 0. The function f(z) = z' 1 ^ - z ) a e- bm T^ z , 
where m\ = m 2 + p\ and m is the mass of the quark, 
expresses the probability that a fraction z of the available 
energy will be carried away by a newly created hadron. 
The parameters (a, b) were modified from their default 
values (0.3, 0.58) to (0.5, 0.1), a change making the frag- 
mentation softer. This modification was inspired by a 
similar study done by the HERMES experiment [23, 
and seems to work also in the present case, with smaller 
deviations from the default values. However, we are look- 
ing here at a particular sample and have no possibility 
to check if the Monte Carlo sample generated with these 
modifications would correctly describe the full data. The 
uncertainty connected with these modifications has been 
estimated and included in the systematic error. 

B. Simulation of experimental conditions 

The scattered muon track of each simulated event was 
followed through the magnet aperture. Trigger condi- 
tions were checked and prescaling factors applied in or- 
der to reproduce the relative trigger rates in the simu- 
lated sample. Kinematic smearing was applied to muon 
and hadron tracks and geometric smearing to the vertex 
position. In addition, the loss of tracks due to chamber 
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inefficiencies was taken into account by applying detector 
plane efficiencies to the simulated events and by removing 
the tracks which did not fulfill the minimal requirements 
for reconstruction. 

Secondary interactions of hadrons have to be taken into 
account to reproduce the distribution of interaction ver- 
tices along the target axis. Hadrons were rejected from 
the sample according to the probability of re-interaction 
in the polarized target material. As an example, Fig. [21 
shows the agreement obtained for the vertex position 
along the beam axis in one of the proton data sets. 

The simulation was performed for each year of data 
taking separately. To get a good description of the kine- 
matic variables it was required to use specific beam pa- 
rameters for every year, including small changes in an- 
gles, and to take into account the exact target position. 



C. Comparison of simulations and data 

The distributions of kinematic variables as well as the 
particle distributions in detectors were checked with iden- 
tical selection criteria applied to data and MC. For the 
simulated events the cuts were applied to the smeared 
variables. The distributions for data and MC were nor- 
malized to the same number of events. The distributions 
of x and Q 2 for interactions on protons are presented in 
Fig. [3] The obtained agreement is at the level of f 0-25% 
for all kinematic event variables. The level of agreement 
for deuterons is very similar p9| . 

The same comparisons were done for hadron variables. 
For simulations performed with the unmodified fragmen- 
tation function clear discrepancies are observed for the 
hadron production angle 9 and the longitudinal momen- 
tum pl, while satisfactory agreement is obtained for pt, 
except at the highest values. The observed differences 
at the highest values of pr can be explained by the ap- 
proximate description of the non-Gaussian tails of the 
distributions used for smearing and by the effects of real 
photon radiation, which are not taken into account in 
the present analysis, ft was checked that the discrep- 
ancy for the 9 angle could not be removed by using dif- 
ferent smearing parameterizations or even by an artifi- 
cial increase of smearing. Agreement between data and 
simulation could only be achieved by applying a cut on 
the hadron production angle 9 > 0.02 rad. This cut, 
however, removes about 25% of the selected sample and 
cannot be justified since there is no reason why the sim- 
ulation should not describe the hadrons produced at low 
9. Therefore modified simulation conditions providing a 
better description of the data were searched for. 

When the modifications of the fragmentation function 
parameters are applied (cf. Section V.A), the agree- 
ment becomes satisfactory over a wide range of 9 and 
Pl- The comparison of the pi, and 9 distributions is 
shown in Fig. 0] for the hadron with highest px- The 
second hadron is also well described by the MC [2|j . We 
concluded that the parameters of the longitudinal frag- 



mentation function f(z) have to be modified in order to 
obtain a good description of the data over the full range 
of hadron production angle 9. Since it is difficult to check 
if the modified set of parameters correctly describes the 
semi-inclusive hadron distributions, the analysis has been 
performed in parallel with modified fragmentation as well 
as with the standard fragmentation and an additional cut 
on 9 > 0.02 rad. 



VI. SELECTION OF THE PGF PROCESS 

In order to compare the merits of various selections of 
PGF events, we will use the efficiency e, which is the 
ratio of the number of PGF events accepted by the se- 
lection criteria to the total number of PGF events, and 
the purity Rpgf (Eq.3), which is the ratio of the number 
of selected PGF events to the total number of selected 
events. The optimal selection is obviously the one pro- 
viding the highest values of e and Rpgf but, in general, 
an increase of the former will result in a decrease of the 
latter. 

The purity is 0.11 for the full sample of events with 
at least 2 charged hadrons. The additional requirement 
of two hadrons with pp > 0.7 GeV defines our reference 
sample for which Rpgf = 0.24 and, by definition, e = 1. 

The effects of cuts were studied for the following vari- 
ables: pti, the sum p\ x + p|, 2 , hadron charges (same 
or opposite sign), the azimuthal angle (j> between the mo- 
menta of the two hadrons with respect to the virtual pho- 
ton direction, and the invariant mass of the two hadrons 
(see also Ref. |3(| ). It was found that the selection on 
J^Pt i s optimal for enhancing the PGF purity and that 
further requirements on the hadron charges do not bring 
any significant improvement. Fig. [3] shows the variation 
of Rpgf with e when the cut on IS varied up to 

4 GeV 2 . It is seen that the purity increases only very 
slowly when the cut is made more restrictive while the 
efficiency drops very rapidly. This can be understood by 
the fact that one of the background processes (QCD-C) 
has a similar dependence on the J^Pt cu * as PGF. The 
approximation made in Eq. @ by the use of A\ for the 
asymmetry on quarks is only valid if the fraction of PGF 
events in the selected sample is much higher than in the 
inclusive one, i.e. close to the maximum value of 0.33. 
The efficiency also needs to be sufficiently high to allow 
a meaningful analysis. As a compromise, we have fixed 
the cut at 2.5 GeV 2 , which corresponds to e = 0.30 and 
Rpgf — 0.31. 

The combination of several variables into a single pa- 
rameter has also been investigated in a classification pro- 
cedure based on a neural network [2^, |^| . We considered 
the variables which characterize the DIS event (x, Q 2 , y, 
and the multiplicity of tracks) and those which describe 
the two selected hadrons with highest pp (transverse and 
longitudinal hadron momenta, charges of the hadrons, 
energy fraction of the hadrons, and the azimuthal an- 
gle <fi). The classification procedure was trained on a 
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Monte Carlo sample where the actual process is known 
for each event. As a result, the procedure provides a sin- 
gle value, called "NN response", within the range (0,1). 
High values of this response correspond to events which, 
according to the classification algorithm, are more likely 
to be PGF than background processes. A threshold on 
the network response can thus be used to select a PGF 
enriched sample. 

The variation of Rpgf vs. e for various choices of 
the NN response threshold is also shown in Fig. [3] It 
is observed that at equal efficiency the NN approach al- 
ways provides samples with higher purity than the selec- 
tion based on 52Pr- For further analysis, a threshold of 
0.26 was chosen, which corresponds to Rpgf — 0.33 and 
e = 0.56. A similar purity is obtained with the 52 Pt cu ^ 
at 2.5 GeV 2 but with an efficiency of 30%. Therefore a 
better statistical precision on the measured asymmetry 
will be obtained with the neural network method. Alter- 
natively, a higher NN threshold corresponding to a PGF 
efficiency of 30% would yield a sample where the purity 
is about 37%, i.e. 6% higher than the value obtained 
with the 52 Pt cu ^- The comparison of the two selected 
samples shows that the NN procedure selects a large frac- 
tion of events with 52 Pt > 2. 5 GeV 2 but also covers the 
lower range of 52 Pt- ^ was & ^ so checked that the dis- 
tributions of NN responses are compatible for data and 
Monte-Carlo events. 



VII. SPIN ASYMMETRIES A eN ~* ehhx 

The SMC data taken from 1993 to 1996 were split into 
periods of data taking, corresponding to about 15 days 
each. The asymmetry for a given year is the weighted 
average of the asymmetries calculated for each period of 
data taking. Splitting the data into smaller subsamples 
gives identical results within the expected statistical fluc- 
tuations. The distribution of the vertex position along 
the beam axis, as presented in Fig. [3 shows that the ra- 
tio of acceptances for the upstream to downstream target 
cells is about 0.7. The method used for asymmetry calcu- 
lation, described in (l8| . is suited for such an acceptance 
difference. 

The asymmetry calculations were done for the entire 
sample which has a purity Rpgf = 0.24 and for the two 
selection methods with enhanced Rpgf (52 Pt > 2.5 
GeV 2 and NN response > 0.26). The results given in 
Fig. and Table |U show that the asymmetries do not 
change significantly with the selection. Also the asym- 
metries obtained for proton and deuteron are compatible 
within errors. The statistical error is larger for the selec- 
tion based on 52 Pt because a smaller fraction of events 
is selected (28 % vs. 42 %). 

The errors of the measured A £N ~* ehhx asymmetry 
for the selected samples are dominated by statistics. 
The contributions to the systematic uncertainty on 
A l N^ihhx are detailed in Table |n| for the two selec- 
tions with enhanced Rpgf- The most significant ones 



come from the false asymmetries, the fraction of radia- 
tive processes (p) and the polarized radiative corrections. 
For the false asymmetries an upper limit from the time 
variation of the acceptance was taken under the assump- 
tion that the reconstruction for each of the three tracks 
(scattered muon and two hadrons) is affected indepen- 
dently. The method used for estimating these effects is 
described in Ref . |30| . The radiative corrections are small 
due to the limited phase space for real photon emission in 
events where a significant fraction of the available energy 
is taken by the two hadrons with large pt- The uncertain- 
ties in p and polarized radiative corrections were taken 
equal to the full size of the inelastic contribution. The 
effect of real photon radiation on the event kinematics 
and, in particular, on the value of pt itself has not been 
taken into account in view of the limited precision of the 
present analysis. 

VIII. DETERMINATION OF THE GLUON 
POLARIZATION 



The gluon polarization is evaluated from Eq. @ us- 
ing the measured A eN ^ ihhx asymmetry, obtained for the 
samples with enhanced R pg f , quoted in Table 1 . In view 
of the strong dependence of the resulting gluon polariza- 
tion on the information obtained from the Monte Carlo, 
special attention was given to the agreement of data and 
simulated events (Figs. 2-4). 

The asymmetry A\{x) for each event is taken from a 
fit to all experimental data and averaged for the full pro- 
ton and deuteron samples. The partonic asymmetries 
&ll for each sub-process are calculated for each Monte- 
Carlo event and averaged. Their averaged values for 
LP and QCD-C are very similar for the two selections 
namely, (a LL ) LP = 0.8 and (a LL )^ CD - c = 0.6. The 
values for PGF are (a LL ) PGF = -0.44 and -0.49 for 
the Epip cut and the NN selection, respectively. After 
selection on Y^p\ the final proton sample consists of 26% 
LP, 43% QCD-C and 31% PGF, while for the neural net- 
work the fractions are Rlp = 38%, Rqcd-c = 29% and 
Rpgf — 33%. The contributions of different processes 
for the proton and deuteron samples differ by less than 
2%. 

The gluon polarization is determined for the kinematic 
region covered by the selected sample and corresponds to 
a given fraction of nucleon momentum carried by gluons 
77: 

v = ^ + i). ( 5 ) 

This quantity is known for simulated events but cannot 
be directly determined from the data. Nevertheless, s 
can be approximately calculated from the virtual photon 
energy in the laboratory system and from the angles (6± , 
62 ) defined by the directions of the two hadrons with 
respect to the virtual photon: 

s « ^ 2 tg0itg0 2 - (6) 



7 



To check the validity of this approximation in our kine- 
matic conditions, we have compared the generated r\ and 
the one calculated from the above equation for selected 
PGF events. The calculated values are on average 25% 
higher than the generated ones. The averaged value of 
the generated r\ for the selected PGF events in the Monte 
Carlo is used as the reference value for the result on 
AG/G. We have also checked the average values of 77 
calculated for all simulated events and obtained the val- 
ues 0.15 for the cut J^Pt > 2 - 5 Gey2 and °- 10 for the 
NN response > 0.26. For both selection methods the val- 
ues of r\ calculated for all simulated events and for data 
are very close. The results on the gluon polarization and 
the values of (rf) are presented in Table ITTll 

In addition to the systematic errors on the measured 
asymmetry discussed in Section 7 and given in Table ITTI 
the asymmetry A\, the fractions R, and the partonic 
asymmetries {(III) contribute to the systematic error on 
AG/G. The contribution due to the asymmetry A\ is 
determined from the uncertainty on A\ at the averaged 
value of x and thus from the errors on the fit parameters. 
The value of A\ at the average x agrees with the average 
A\ calculated from the fit for each event to within 0.001. 

The dominant contributions to the systematic error are 
due to the uncertainties on the values of R and (clll)- 
They are estimated by comparing the results obtained 
from Monte Carlo simulations with different parameters. 
For this purpose, a sample of LEPTO events was gener- 
ated with the same kinematic and hadron selections but 
with modified renormalization and factorization scales, 
cut-offs and fragmentation function parameters. Scales of 
Q 2 /2 and 2 Q 2 were used for comparison and provide an 
estimate of the stability of the leading order approxima- 
tion used here. Results with standard and modified pa- 
rameters (see Section 5.1) in the fragmentation function 
were compared. Since only the simulations which repro- 
duce the data should be considered, a cut on the hadron 
angle 9 was applied, as explained in Section 5.3. The 
value of the gluon polarization calculated with this new 
Monte Carlo sample was compared to the one obtained 
under the conditions described in Section 5.1. This pro- 
cedure was repeated several times with slightly different 
cuts and with different neural network thresholds. For 
the neural network the procedure is complicated by the 
fact that any change in the simulation procedure leads to 
a different selection on the data. To avoid the fluctuation 
of the gluon polarization due to variation of the measured 
asymmetry, the value of this asymmetry was artificially 
frozen when comparing results for different MC samples. 
The individual contributions to the systematic error are 
given, for both selection methods, in Table IIVI It was 
checked that the effect of combined modifications in the 
Monte Carlo is smaller than the sum of the individual un- 
certainties. The maximal variation of Rpgf and {olll) 
was found to be 20% and 4% respectively. 

As discussed before, the neural network selection pro- 
vides a more accurate result than the selection based on 
Hp\ cuts. However, the statistical error is too large to 



draw definitive conclusions on the contribution of AG to 
the nucleon spin. The systematic uncertainty is small 
compared to the statistical error. The demand of a good 
agreement of the simulation with the data sets an impor- 
tant limitation on the estimated systematic uncertain- 
ties. For this reason, an increase in statistical precision 
is expected also to lead to further improved systematic 
uncertainty estimates. 

Averaging the results for proton and deuteron obtained 
with the neural network classification we obtain AG/G = 
-0.20±0.28±0.10. 



IX. CONCLUSIONS 

We have evaluated for the first time the gluon polar- 
ization from the spin asymmetries measured in lepton- 
nucleon DIS events with Q 2 > 1 GeV 2 including two 
hadrons with large transverse momentum in the final 
state. The analysis is performed at leading order in QCD 
and based on the comparison of selected data samples 
with simulated events provided by the LEPTO genera- 
tor. The partonic asymmetry clll is mostly negative for 
the photon-gluon fusion process while it is positive for 
the two competing processes, leading process and gluon 
radiation. The relative contribution of photon-gluon fu- 
sion is enhanced to about 30% by applying a cut on 
Ep|, > 2.5 GeV 2 or by using a neural network classi- 
fication. 

The average gluon polarization obtained for the SMC 
data is close to zero with a large statistical error (~ 0.30). 
The accuracy is limited by the reduction to less than 1% 
of the DIS sample by the hadron selection requirements. 
It is thus expected to be improved by higher counting 
rates and larger hadron acceptance in ongoing and future 
experiments. 
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b) 




y* 



FIG. 1: Lowest order diagrams for DIS 7* absorption: a) leading process (LP), b) gluon radiation (QCD-C), c) photon-gluon 
fusion (PGF). 
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FIG. 2: Distribution of vertices along the beam axis. Points correspond to the proton data from 1993 and the histogram to 
the corresponding MC simulation. 
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FIG. 3: The x and Q 2 distributions for the proton case: points correspond to the data and histograms to the Monte Carlo 
simulation. 
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FIG. 4: Distributions of longitudinal momentum and scattering angle for the hadron with the highest pr- Points correspond 
to the proton data collected in 1993, histograms to the Monte Carlo simulations with the modified fragmentation function. 
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FIG. 6: Measured asymmetry A ~* , for proton and deuteron, for events with pti,2 > 0.7 GeV cut and after additional 
selections on J^Pt and neural network threshold to increase the purity. 
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TABLE I: Measured cross-section asymmetries A eN ^ ehhx f or proton and deuteron events with pri,2 > 0.7 GeV and in the 
samples selected with the ^2 Pt cu t and with the neural network response threshold, each given with statistical and systematic 



Selection 



AtN-tlhhX 



A iN-*ihhX 



All 

J^Pt > 2 -5 GeV 2 
NN response > 0.26 



0.041±0.037±0.011 
0.018±0.071±0.010 
0.030±0.057±0.010 



0.063±0.050±0.011 
0.054±0.093±0.008 
0.070±0.076±0.010 



TABLE II: The contributions to the systematic error of A lN ^ lhhx with the ^Pt >2.5 GeV 2 cut and with the neural network 
response >0.26 for SMC proton and deuteron data. The first and last contributions are additive; the others are proportional 
to the asymmetry. 



Contributions to the 
systematic error on A m ~* lhhx 



Y,p 



proton data 



NN 



deuteron data 



NN 



False asymmetries 
Target polarization 
Beam polarization 
Dilution factor 

Target composition 

p factor 
Polarized rad. corr. 



0.0049 
0.0005 
0.0007 

0.0003 
0.0018 
0.0083 



0.0049 
0.0008 
0.0011 

0.0001 
0.0030 
0.0083 



0.0044 
0.0016 
0.0021 

0.0002 
0.0054 
0.0020 



0.0044 
0.0023 
0.0029 

0.0001 
0.0076 
0.0020 



Total systematic error 



0.0098 



0.0102 



0.0077 



0.0097 



TABLE III: Gluon polarization for proton and deuteron for the Spy cut and the neural network selection. 



Selection 



< a ), 



' AG \ 
. G ) d 



(V) 



J2pt> 2.5 GeV 2 
NN response> 0.26 



0.11±0.51±0.12 
-0.06±0.35±0.10 



-0.37±0.66±0.12 
-0.47±0.49±0.10 



0.09 
0.07 
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TABLE IV: Contributions to the systematic error on gluon polarization for two methods of event selection. 



Source of the uncertainty Ep^ NN 



systematic error 

on A tN ^ thhx 0.072(p) 0.057(d) 0.061(p) 0.063(d) 

precision of A 1 fit 0.042(p) 0.042(d) 0.026(p) 0.028(d) 
scale change 

from Q 2 /2 to 2 Q 2 0.008 0.010 

fragmentation paramr. 0.036 0.034 

cut-offs in matrix elem. 0.015 0.008 



